Entanglement evolution of a two-qubit system with decay beyond rotating-wave 

approximation 
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Two noninteracting atoms, initially entangled in Bell states, are coupled to a one-mode cavity. 
Based on the reduced non-perturbative quantum master equation, the entanglement evolution of 
the two atoms with decay is investigated beyond rotating-wave approximation. It is shown that the 
counter-rotating wave terms have great influence on the disentanglement behavior. The phenomenon 
of entanglement sudden death and entanglement sudden birth will occur. 
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I. INTRODUCTION 

Entanglement plays a central role in the field of quan- 
tum information science where it is considered as a valu- 
able resource for some non-classical tasks, such as: quan- 
tum computation]!], 0, S| , quantum teleportation [3, su- 
perdense codingU and quantum cryptography [g, Q. 
Applications of interest have triggered research on the 
dynamical behavior of entanglement in order to control 
quantum disentanglement [ct M, [Tpj . 

Recently, the dynamics of entanglement in bi-partite 
systems has been under extensive research [IJ [12, lL^.ll4l . 
13 EE [3 El, [3 III since the work of Yu and Eberly|2lj], 
in which the entanglement may terminate abruptly in a 
finite time while coherence is lost asymptotically. This 
effect is termed entanglement sudden death(ESD). Sub- 
sequently, C.E.Lo'pez etc. put forward a new term "en- 
tanglement sudden birth" (ESB) [22] and try to present 
an explanatory study of multipartite entanglement evo- 
lution. 

In the previous studies, the rotating-wave approxima- 
tion(RWA), which neglecting counter rotating terms cor- 
responding to the emission and absorption of virtual pho- 
ton without energy conservation, is widely used. Gener- 
ally, the coupling ratio of the atom-field is of the order 
1CP 7 ~ 1CP 6 in atom-field cavity systems and the RWA 
is justified. Recently, D.Meiser etc. have investigated the 
cavity systems with superstrong coupling 23]. It can be 
seen that the effect of counter-rotating terms should be 
considered in the strong coupling regime, such as solid 
state systems. 

Different from the previous works, we study the disen- 
tanglement between a pair of qubits with decay which 
interacting with a cavity beyond the RWA. Consider- 
ing the real processing, the atomic decay is unavoidable, 
which can be caused by two different physical mecha- 
nisms. The first is the transition without light radiation. 
In this process the energy is emitted by thermal energy 



or the other forms. The second is the radiation transi- 
tion. There will emit photons with the atomic transition. 
In our model the second is the main physical mechanism 
because the collision probability between the two atoms 
is very small. Thus the discussion upon the effect of 
atomic decay caused by spontaneous emission to the en- 
tanglement evolution is necessary. 

The paper is organized as follows. In section II, the 
two-qubit model with decay is presented and the reduced 
non-perturbative quantum master equation of atoms is 
derived. In section III, the entanglement evolution of two 
initially entangled atoms is investigated by using Woot- 
ters' concurrence [24J. In the last section, the conclusions 
are given. 



II. MODEL 

Consider two noninteracting atoms A and B, which 
are interacting with a single-mode cavity resonantly. For 
simplicity, we assume that the two atoms have the same 
parameters. It can be described by the following Hamil- 
tonian: 



where 



H = H a +H f + H af , 
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Hr = coa a, 



H a f = Xa(a + + a), 
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(4) 



where a = £f =1 (cr+ + a i ), a z% = |e,-)(ei| - 
err = |<7j)(ej| are the atomic operators with 
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\gi)(gi\ and 
e;) and 

being the excited and ground states of the ith atom; lu is 
the atomic transition frequency between the ground state 
and the excited state; A is the coupling constant between 
atom and cavity; a and denote the annihilation and 
creation operators of the cavity field mode corresponding 
frequency w. 
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Taking the atomic spontaneous emission into consid- 
eration, the reduced non-perturbative quantum mas- 
ter equation of atoms could be obtained by path 
integrals [25j. 



^p a = - iCaPa- J Q ds(C af e- lC ^-^ Cafe-^-ty 

where Cq, C a and C a f are Liouvillian operators defined 
as 



Here we restrict our analysis to the initial entangled 
states 

\4>) = or|01) + 0\1O), |V) = a|00) + /5| 11) (11) 

where a is real, /3 = \/3\e lS and a 2 + \/3\ 2 = 1. For these 
fP a two entangled states, the initial atomic density matrix 
has an "X" structure [26[ which is maintained during the 
(5) evolution (27}. The reduced density matrix of the two 
atoms p a , in the standard product basis B = {|1) = 
|ee> , |2> = \eg),\3) = |se),|4) = \gg)}, could be written 
as 



£qP = [H a + Hf,p], 
C a p = [H a ,p], 
C af p=[H af ,p], (6) 

and < • • • > f stands for partial trace of cavity mode; 7 
is the atomic decay constant. 

Assuming that the cavity field is initially in vacuum 
state, the non-perturbative reduced master equation of 
the atoms could be derived from Eq.(5), 



d_ 

dt 



Pa 



aX 2 a[cr+ +a+,p] - f\ 2 (j[a^ + u^,p] 
a*A 2 [crf + 02 , p](T + /*A 2 [(T ] f + &2 , p](J 
7rf tn _-_ _+ _+ - „ + 



where 



Pa - PaOja i )(7) 

1 — exp(-iAt) 



iA 

exp(iSt) 



1 



iS 

and A = uj + u>o, S = luq — to, /*is conjugate of /. 



III. DISENTANGLEMENT 



A. Initial state and the entanglement measurement 



In this section, we will use Wootters' concurrence to 



quantify the degree of entanglement 24} ■ For two qubits, 
the concurrence is calculated from the density matrix p 
for qubits A and B: 



C(p) = max(0, v/Ai - V\ 2 - V A 3 - V A 4), (9) 
where Xi are the eigenvalues of the matrix 

Q = PAB{tf®Vy)p* AB {°t® u y) ( 10 ) 
arranged in decreasing order. Here p* AB denotes the com- 
plex conjugation of pab, and (Jy^ is the standard Pauli 
matrix acting in the space of qubit A (or B). The con- 
currence varies from C{p) = for an unentangled state 
to C(p) = 1 for a maximally entangled state. 



Pa 
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(12) 



The concurrence of p a could be obtained 



cm = 



max{0, 2\p2 3 \ 
- max{0, 2|p 14 | 



2V/°ll/°44}, 
" 2VP22P33}- 



(13) 



B. Numerical results and discussion 



In order to study the counter rotating terms effect on 
the system, we investigate the entanglement evolution of 
two atoms by numerical calculation. For simplicity, the 
resonant case uj = loq is considered. 

First, we focus on the disentanglement of two qubits 
with the initial state of \cj>}. The entanglement evolution 
of two decayed atoms is shown in Fig.l, Fig. 2 and Fig. 3 
with different parameters. One could see that the change 
(g) of against a 2 is symmetrical because of the symmetry 
of the initial state \<j>). 

When uj = A, Fig.l shows that the concurrence 
changes with the initial value a 2 and the time Xt. It can 
be seen that the concurrence decreases to zero in a finite 
time and the entanglement undergoes the so-called ESD. 
The tendency is similar to that of the above case when 
lu = 3A, as shown in Fig. 2. 

As uj = 10A, Fig. 3 reveals the time evolvement of con- 
currence. The tendency is different from that of the 
above cases. It is easy to note that the entanglement 
decays to zero in a finite time, then revives with small 
amplitude and disappears permanently at last. Partic- 
ularly, when the two atoms are not entangled initially, 
namely a 2 = or a 2 = 1, there is entanglement sudden 
birth and its amplitude is larger than that of the case 
< a 2 < 1, then the entanglement decreases to zero 
eventually due to the atomic decay. 

Alternatively, we focus on the disentanglement of two 
qubits with the initial state \ip). The entanglement 
evolvement is investigated in Fig. 4, Fig. 5 and Fig. 6. The 
figures are not symmetrical to a 2 because the initial state 
of the two atoms is asymmetrical. 

When lu = A, Fig. 4 shows that the concurrence de- 
cays exponentially to zero in almost all cases except that 



FIG. 1: The concurrence Cp as functions of a 2 and the time FIG 4 . The concurrence Cil as functions of a 2 and the time 
At with u = A. 6 i/2 Xt with LU = y 




FIG. 2: The concurrence C*0 as functions of a 2 and the time 
Aa£ with uj = 3A. 



the value of a 2 is near 0. In the case of small a 2 , there 
exists small fluctuation before C^p vanishes permanently. 
Particularly, when a 2 = 0, there exists small ESB and 
ESD. This effect is resulted from the coaction of rotating 
wave process and counter-rotating process on the whole 
system. 

By contrast, as u> — 3A, Fig. 5 reveals that the ampli- 
tude of fluctuation is higher than that of uj — A in the 
case of small a 2 . It is also found that the range existing 
the fluctuation is larger than that of the previous case. 
But there is no ESB in the case of a 2 = 0. 




X t 



FIG. 5: The concurrence Cp as functions of a 2 and the time 
Xt for u) = 3A. 



When uj = 10A, the time evolution of the concurrence 
is plotted in Fig. 6. The concurrence decreases monoton- 
ically and exponentially to zero when a 2 is near 1. The 
smaller the a 2 is, the more intense the fluctuation is. The 
concurrence disappears permanently at last because the 
atoms decay. 




FIG. 3: The concurrence Cj, as functions of a and the time FIG. 6: The concurrence as functions of a 2 and the time 
X A t with u = 10A. \t for u) = 10A. 
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IV. CONCLUSIONS 

In summary, we have investigated the disentanglement 
of two noninteracting atoms with decay coupling to a one- 
mode cavity resonantly beyond RWA. It is shown that 
the entanglement evolution is dependent on the ratio of 
the atom-field coupling divided by the atomic transition 
frequency and the phenomena of entanglement evolution 
are rich. The physical mechanism behind the phenom- 
ena is the process of emission and absorption of virtual 
photon and the atomic decay. 

The study of entanglement evolution beyond RWA is 
a significant problem because of its importance to the 
field of strong coupling. It will help ones to deal with the 
practical case of solid state system. The consideration of 
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